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The effective string theory emerging from the bilocal approximation to the Method of Vacuum 
| Correlators in gluodynamics is shown to be well described by the 4D theory of the massive Abelian 
Kalb-Ramond field interacting with the string, which is known to be the low-energy limit of 
the Universal Confining String Theory. This correspondence follows from the agreement of the 
behaviour of the coefficient functions, which parametrize the gauge-invariant correlator of two 
gluonic field strength tensors, known from the lattice data, with their values obtained from the 
propagator of the Kalb-Ramond field. We discuss this correspondence in several aspects and 
demonstrate that the mass of the Kalb-Ramond field in this approach plays the role of the inverse 
5_i ■ correlation length of the vacuum, so that in the massless limit string picture disappears. Next, we 



apply the background field method, known in the theory of nonlinear sigma models, to obtain the 
action, which is quadratic in quantum fluctuations around a given (e.g. minimal) string world- 
sheet. Several nontrivial types of couplings of these fluctuations with the background world-sheet 
are obtained and discussed. 



1. Introduction 

Recently, a new approach to the string representation of the confining phase of gauge theories 
was proposed 1 . This is the so-called Universal Confining String Theory (UCST), which is the 
theory of an Abelian antisymmetric tensor field with a nonlinear action, interacting with the 
string. This field effectively substitutes an infinite number of monopoles in the 3D compact 
QED. It was argued in Ref. 1 that the summation over the branches of the UCST action should 
correspond to the summation over the string world-sheets. It was also proved in Ref. 1 that the 
UCST partition function, which is nothing else but the Wilson average in the 3D compact QED, 
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satisfies loop equations 2 modulo contact terms. However, the statement made in Ref. 1 concerning 
the universality of the wave operator standing on the L.H.S. of the usual loop equations obtained 
from the Yang-Mills theory, which were derived and investigated in Ref. 2, was only a conjecture, 
and the relation of these equations to the equation of motion of the tensor field in the UCST 
was absolutely unclear. In order to clarify this conjecture, in Ref. 3 the loop equation for the 
4D UCST partition function was derived and investigated. In particular, it was demonstrated in 3 
that the wave operator of the obtained loop equation is quite different from the usual one. Also 
the corresponding contact term was calculated explicitly. 

The low-energy limit of the UCST, in which the Wilson loop could be evaluated exactly, was 
discussed in Ref. 1 and investigated in the 4D case in Ref. 4, where the string tension of the 
Nambu-Goto term and the coupling constant of the rigidity term were calculated. The latter 
one occured to be negative, which means that the obtained (Euclidean) string effective action 
is stable 5 . Also the 4D UCST action was derived in Ref. 4 by performing the exact duality 
transformation, while in Ref. 1 it was done only at the semiclassical level. 

However, as it was first already pointed out in Ref. 6, the non-Abelian generalization of the 
above described antisymmetric tensor theory is quite difficult. 

Another, more phenomenological, but on the other hand adapted to the non-Abelian case 
approach to the problem of the string representation of the confining phase of gauge theories was 
developed in Ref. 7-9. It is based on the Method of Vacuum Correlators 10 (MVC) and allows one 
to get the information about the gluodynamics string effective action from the expansion of the 
Wilson average, considered as a statistical weight of this string theory, in powers of the correlation 
length of the vacuum. In this way, in Ref. 7 the sign of the rigidity term was found to be negative 
according to the present lattice data, which means that the MVC is in agreement with the dual 
superconductor model of confinement due to Ref. 5. In Ref. 8, the string effective action obtained 
in Ref. 7, was applied to the derivation of the correction to the Hamiltonian of the QCD string with 
quarks 11 due to the rigidity term. However, up to now it is not clear how the summation over the 
world-sheets could appear within this approach. The most appropriate way in this direction might 
lie in the accounting for the perturbative gluons, which propagate inside the Wilson loop and, as 
it was argued in Ref. 12, should generate world-sheet excitations. In Ref. 9, this conjecture was 
elaborated out by virtue of the integration over perturbative gluons in the Wilson average, which 
in the lowest order of perturbation theory leads to the interaction of the world-sheet elements via 
the exchanges of perturbative gluons, propagating in the nonperturbative background. Expanding 
this interaction in powers of the derivatives w.r.t. the world-sheet coordinates, we finally get 
in the lowest order of this so-called curvature expansion some definite correction to the rigid 
string coupling constant, while the string tension of the Nambu-Goto term acquires no additional 
contribution due to perturbative gluons and keeps its pure nonperturbative value. 

In this Letter, we shall propose an alternative method of description of fluctuations of the 
gluodynamics string world-sheet. To this end, we shall first demonstrate in the next Section that 
the gluodynamics string emerging within the bilocal approximation to the MVC could be effec- 
tively described by the same action of the massive Kalb-Ramond field interacting with the string, 
which describes the low-energy limit of the UCST. This observation will be proved by the cal- 
culation of the coefficient functions standing at two Kronecker structures (which describe surface 
and boundary terms in the string effective action) in the propagator of the Kalb-Ramond field, 
and further comparison of them with the values of the coefficient functions standing at the cor- 
responding structures in the gauge-invariant correlator of two gluonic field strength tensors. The 
latter ones are unfortunately yet not found from the gluodynamics Lagrangian, while equations 
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for correlators were obtained in Refs. 13 and 14 and then investigated in Ref. 15 by making use of 
the stochastic quantization method (see also Ref. 16, where alternative equations for correlators 
following from the non-Abelian Bianchi identities, which were proposed in Ref. 17 and generalized 
in Ref. 13, were investigated). That is why, we shall compare the coefficient functions, standing in 
the propagator of the Kalb-Ramond field, with the lattice data 18 concerning the behaviour of the 
corresponding coefficient functions in gluodynamics (see also Ref. 19, where these functions were 
measured in QCD with dynamical fermions). We shall see that the mass of the Kalb-Ramond 
field plays the role of the inverse correlation length of the vacuum, so that in the sum rules' limit, 
when the correlation length of the vacuum tends to infinity (see the last Ref. in 10), which in our 
model corresponds to the case of the massless Kalb-Ramond field, the string picture disappears, 
and we are left with the boundary terms only. The points described above will be the topics of 
the next Section. 

In Section 3, we shall use the model of the gluodynamics string proposed in Section 2, to 
describe string world-sheet excitations. Correspondingly, for the low-energy limit of the UCST this 
will be an exact procedure rather than a model dependent approach. The world-sheet excitations 
will be described with the help of the background field method developed in Ref. 20 for nonlinear 
sigma models. Namely, we shall split the world-sheet coordinate into a background and quantum 
fluctuations, after which upon the integration over the Kalb-Ramond field we shall derive for the 
latter ones a quadratic action, which contains several nontrivial couplings of quantum fluctuations 
with the background world-sheet with- and without derivatives. 

The main results of the Letter are summarized in the Conclusion. 

In the Appendix, we perform rather nontrivial integration over the Kalb-Ramond field in the 
expression for the 4D UCST partition function. 

2. A Unified Description of the Gluodynamics String and the Low- 
Energy Limit of the UCST 

The partition function of the 4D UCST (which is nothing else but the Wilson average in the 
Euclidean 4D compact QED) in the low-energy limit has the form 1 ' 4 

(W(C)) = N J DB^exp 

Here 

H^x = d^B uX + d v B^ + d\B^ v 
is a strength tensor of the field B^, 

T^ix) = J da„ u (x(^))5(x - x(g)) 
is the vorticity tensor current, e is a dimensionless coupling constant, and A = ^f-\/z, where 

const . 

z ~ e e 2 , and A is a cutoff which is necessary in 4D. Gaussian integration in (1) is carried out 
in the Appendix and yields the UCST low-energy action in the form 

Sucst = J da Xv (x) j da^ p (x') (B Xv (x)B w (x')) . (2) 

Here 



/ 



dx\ 



——H 2 

12A 2 ^ A 



(1) 
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[B Xl/ (x)B w (0)) = {BUx)B»M) (1) + {Bxu{x)B^)Y Z) , 



(3) 



where 



{B Xv (x)B w (Q)y = — -j \Sx^S vp - S^dxpJ, 



(4) 
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In Eqs. (4) and (5), m = ^ is the mass of the Kalb-Ramond field, AVs, i = 0, 1,2,3, stand for 
the Macdonald functions, and one can show that the term 

J da Xv {x) j da^x') {B Xv {x)B w {x')) {2) 

could be rewritten as a boundary one, due to which one can immediately establish a correspondence 
between (W(C)) = exp (— Sucst) and the Wilson average written within the bilocal approximation 
to the MVC. This correspondence yields the following values of the coefficient functions D and 
Di, which parametrize the gauge-invariant correlator of two gluonic field strength tensors 



2 2 

m x 



e m Ki(m \x\) 
8n 2 \x\ 



(6) 



and 



e 2 m 
An 2 x 2 



Ki(m\x\) m ( „ , , lx r , , i IX N 
- ' u + — ( K (m\x\) +K 2 (m\x\) 



(7) 



Notice, that the derivative (or curvature) expansion of the action (2), which in this case is 
equivalent to the ^-expansion, was performed in Ref. 4 and up to a common constant positive 
factor has the form 



SucsT = „, A , / ^ - \ j <* w (w + if (£) I < 8 > 

Here g a b and t^ v stand for the induced metric and the extrinsic curvature tensor of the string 
world-sheet respectively, 

f(y) = J jW), 

y 
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Xn(s) in the last term on the R.H.S. of Eq. (8) parametrizes the contour C in Eq. (1), £^(0) = 
and we have omitted the full derivative term of the form / d 2 ^^/gR, where R is a scalar 
curvature of the world-sheet. As it was already mentioned in the Introduction, the coupling 
constant of the rigidity term in this expansion is negative, which means the stability of the string. 

Asymptotic behaviours of the functions D and D 1 following from Eqs. (6) and (7) at |x| -C ^ 
and \x\ ^> — read 

^ e 2 m 2 . , 

Dl - i^W (10) 



and 



e 2 m 4 e m ^ 



D^^—^ s-, (11) 

8V27T2 (jn \x\) 2 

e 2 m 4 Q~ m \ x \ 

Dl ~ > A Hi*, , ^ (12) 
4V27T2 [m \x\) 2 

respectively. 

Let us comment on the asymptotic behaviours (9)-(12). Firstly, one can see that Eq. (10) is 
in agreement with the MVC, where due to Ref. 10 at the distances which are much smaller than 
the correlation length of the vacuum T g , 

3tv (x z ) 

The difference in the numerical constants in the asymptotic behaviours (10) and (13) is due to the 
colour factor in the one-gluon exchange diagram, which contributes into Eq. (13). This factor, 
which is absent in the asymptotics (10) of the Abelian propagator, could be accounted for by the 
proper tuning of the charge e of the Kalb-Ramond field, if we replace the running a s (x 2 ) in Eq. 
(13) by some fixed value, say approximate a s (x 2 ) by its "frozen" value, which it acquires at the 
confinement scale 21 -4=, where a is the string tension of the Nambu-Goto term. 

However Eq. (9) agrees with MVC only in the lowest order of perturbation theory, when the 
perturbative part of the function D is absent, and D tends to the gluonic condensate at \x\ — > 0, 
i.e. this agreement is only in a sense that D <^ D-y. One can see that Eq. (9) does not coincide 
with the behaviour of the function D (x 2 ) of the type 



* (aln(W)+A 



where a, (3 stand for some constants, and M is a certain mass parameter, which takes place in 
the next-to- leading order of perturbation theory 22 . This is due to the fact that such a behaviour 
is a specific property of non- Abelian theories 23 and presumably could not be reproduced by any 
Abelian model. 

Secondly, comparing asymptotic behaviours (11) and (12) we see that the function D\ falls off 
much faster than the function D, which is in agreement with the lattice data 18 . The exponential 
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falls-off of the functions D and D\ also agree with Ref. 18, while the preexponential power-like 
behaviours do not. However the 7-Lj-fLts of these behaviours used in Refs. 18 and 19 were 

motivated by the short-distance asymptotics of the function Di, which according to Eq. (10) is 
reproduced in our model as well. 

Eqs. (11) and (12) tell us once more that the mass m of the Kalb-Ramond field in our approach 
should be treated as an inverse correlation length of the vacuum T g , so that in the "string limit 
of QCD" 24 , when the string tension of the Nambu-Goto term, a, is kept fixed at vanishing T g , we 



have a correspondence of the type m ~ y m the opposite regime when the Kalb-Ramond field 
is massless (or equivalently, in the strong coupling regime of the theory (1)), which corresponds 
to the QCD sum rules' case (see the last Ref. in 10), the "confining" part (4) of the propagator 
of the Kalb-Ramond field vanishes, and we are left only with the "boundary" part (5), whose 
contribution to Sucst takes the form 



J da Xv {x) j da w {x') (B Xv (x)B fip (x')) {2) = ^ j dx^ j dx'^ 1 , 

n c c ^ ' 

which could be anticipated from the very beginning. This result is also in agreement with Ref. 
4, where it was shown that in the strong coupling regime of the theory (1) strings are completely 
suppressed when A — > +00. 

3. Description of the World-Sheet Excitations 

In order to describe fluctuations of strings in the model (1), we shall adopt the background- 
field method developed in Ref. 20 for the nonlinear sigma models. The three differences of our 
case with the one considered in Ref. 20 are the absence of the Polyakov term in the action, 
flatness of the field manifold (i.e. the space-time), since we consider an arbitrary Wilson loop, not 
necessarily lieing on the unit sphere, and the necessity of the eventual integration over the field 
B^. Consequently, the geodesies passing through the background manifold y M (£) and the excited 
manifold £ M (£) = + where stands for the world-sheet fluctuation, are straight, 

Pm(£' s ) = Vfj.i.0 + S2; m(0> where s denotes the arc-length parameter, < s < 1. The expansion 
of the action (2) in powers of quantum fluctuations z M (£) could be performed by virtue of the 
following generating functional, which is the arc-dependent term describing the interaction of the 
Kalb-Ramond field with the string in Eq. (1), 

/ [p(£, s)]=if d\B, u [p(£, s)\ e ab (d aP ,(t s)) (d bP ^, s)) . (14) 
Then the term containing n quantum fluctuations reads as 

/(»)= a )] , 

s=0 

and we get from Eq. (14) 

J<°> = i j da, u (y(0)B,u [y(01 , (15) 

/W = i J da, u {y{())z x {i)H, uX [y(0] , (16) 
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and 



1^ = i J rfWOe" 6 (d a y,(0) ((d b z x (0) H v , x [y(0] + \z a (0 (d b y x (0) d a H u , x [y(0]) , (17) 



where during the derivation of Eqs. (16) and (17) we have omitted several full derivative terms. 

It is worth mentioning, that as it was discussed in Ref. 20, the terms (15)-(17) are necessary 
and sufficient to determine all one-loop ultraviolet divergences in the theory (14) at s = 0. This 
statement is important for the model of the gluodynamics string with the action (2), proposed 
in the previous Section, since in the UCST case, Eq. (2) describes only its low-energy effective 
action, which does not suffer from the ultraviolet divergences. It is also important for the Abelian 
Higgs Model in the Londons' limit 25 , where one can show 26 that the partition function (1) with 
an additional integration over metrics and T^ u (x) corresponding to a closed surface is nothing else 
but a 't Hooft loop average defined on the string world-sheet. 

In order to get the desirable action, quadratic in quantum fluctuations z M (£)> we shall first 
carry out the integral 



JV I £>£!„„ exp - J dxi-^H] 



It occurs to be equal to 



exp 



- / da Xu (y(0) I da^ivH')) \{B Xv [y{i)]B w + 2z a (£)- 







dy*(0 



(B x Ay(0}B»Ma) {1) + Za(Oz^> 



d 2 



77 {B\v m)\B„w)\) 



(18) 



where from now on we shall omit all the boundary terms. 

Secondly, one should substitute the saddle-point of the integral 

into Eq. (17). This saddle-point reads 



B^- [y(0] 



■ 9 ^ 

ie m 
2?r 2 



\y(()-y((')\ 

and upon its substitution into Eq. (17), accounting for Eq. (18), and making use of Eq. (4), we 
finally get the following value of the action quadratic in quantum fluctuations 



S, 



quadr. 



2 S 

em 



/ *v (v(fl) J 



+m (K + K 2 )\ + 



\y(0 - v(0 



*«(fl(y(fl-y(Q)< 
lv(fl-v(OI 

r ( K l 



2K X 
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|y(0 - y(0 
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m,„ ...\ (»(«-!/«'))«(!/(«-!/«'))«, (Jm> 



__ (A - + A , ) + _ /<3 



where 



^ ( / ^,(0 V* [y(0] e ab (d a y,(0) + 

+ ^ / *V (2/(0) ^(0^(0^a [2/(0]) > (19) 
[2/(0] = 



J da, x (y(0) (z/(0 - 2/(0)1 + / (y(0) (?/(0 - 2/(0)) A + / (?/(0)) (2/(0 - 2/(0)), 



(2/(0 -2/(0)Hl2/(0- 2/(0)1 " 2 - " ' (20) 

and everywhere in Eqs. (19) and (20) the arguments of the Macdonald functions are the same, 
m 1 2/(0 — 2/(0)1- We have also kept the pure background term on the R.H.S. of Eq. (19). The 
other terms yield the couplings of quantum fluctuations with the background world-sheet. All of 
these terms except one do not contain the derivatives of quantum fluctuations. 

4. Conclusion 

In this Letter, we have modelled the gluodynamics string effective action by the theory of the 
massive Abelian Kalb-Ramond field, interacting with the string. The partition function (1) of this 
theory is nothing else but the low-energy expression for the partition function of the UCST, which 
is in fact the low-energy limit of the Wilson average in the 4D compact QED, rewritten in terms 
of the dual antisymmetric tensor field. This observation provides us with a unified description of 
the gluodynamics string and the UCST. 

An approach to the gluodynamics string with the help of the theory (1) has been justified 
in Section 2, where it has been shown that the small- and large-distance asymptotic behaviours 
of the coefficient functions, which parametrize the gauge-invariant correlator of two gluonic field 
strength tensors within the MVC and could be extracted from the lattice measurements, are in a 
good agreement with the ones of the corresponding functions standing in the propagator of the 
Kalb-Ramond field. It has also been demonstrated that the mass of the Kalb-Ramond field plays 
the role of the inverse correlation length of the vacuum, so that in the massless limit strings are 
suppressed, which is in agreement with Refs. 4 and 24. 

In Section 3, we have applied the background field formalism, developed in Ref. 20 for the 
nonlinear sigma models, in order to describe fluctuations of strings in our model. To this end, 
we have splitted the string world-sheet coordinate into a background part, corresponding to some 
given world-sheet (e.g. one with the minimal area), and quantum fluctuations, and using the arc- 
dependent string action as a generating functional, performed the expansion of the term, which 
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describes the interaction of the string with the Kalb-Ramond field, up to the second order in 
quantum fluctuations. Finally, upon the integration over the Kalb-Ramond field, we have derived 
an action, which is quadratic in quantum fluctuations and contains the pure background part and 
the terms describing the interaction of the background world-sheet with the quantum fluctuations. 
This action is given by formulae (19) and (20) and describes fluctuations of strings in the UCST 
and in our model of the gluodynamics string. 

In the forthcoming paper 26 , we shall demonstrate that there in fact exists a physical connection 
between UCST and gluodynamics based on the Abelian Higgs Model. 
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Appendix. Integration over the Kalb-Ramond Field in the UCST 
Partition Function (1) 

According to the general rule, in order to calculate the Gaussian integral (1) one should substi- 
tute its saddle-point value back into the integrand. The saddle-point equation in the momentum 
representation reads 

This equation can be most easily solved by rewriting it in the following way 

(p 2 Px», a p + m 2 l Xu ^) B$-{p) = 2zA 2 T Xu (p), (A.l) 

where we have introduced the following projection operators (see, for example, Appendix to the 
last Ref. in 15 ; functional generalization of these operators has been introduced in Ref. 3) 



and 



1 



where = 8^ — ^ iL . These projection operators possess the following properties 

^-pu,Xp ^-up,Xp ~^-pv,p\ ^-Xp,pvi (^'2) 

^-p,v,\p^-\p,aj3 = lpu,a/3 (^-3) 

(the same properties hold for P^Ap), and 

P,,,a p (1 - P)^ = 0. (A4) 
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By virtue of properties (A.2)-(A.4), the solution of Eq. (A.l) reads 



2iA 2 



p2 _|_ m 2 



1 + ^(1-P) 

m 2 



Xu,a/3 



which, once being substituted back into partition function (1), yields for it the following expression 



(W(C)) = exp 



-A* 



dp 



(27r)V + m 2 



1 + ^(1 

m 2 



T^(-p)T af3 (p) 



(A.5) 



Rewriting Eq. (A. 5) in the coordinate representation we arrive at Eq. (2). From Eq. (A. 5) one 
can easily see that the term on its R.H.S. proportional to the projection operator 1 — P yields in 
the coordinate representation the boundary term, i.e. Eq. (5), whereas the term proportional to 
the projection operator 1 yields Eq. (4). 
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